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Abstract 

The role of various magnetic inter-chain couplings has been investigated recently by numerical methods in doped 
frustrated quantum spin chains. A non-magnetic dopant introduced in a gapped spin chain releases a free spin- 
1/2 soliton. The formation of a local magnetic moment has been analyzed in term of soliton confinement. A 
four-spin coupling which might originate from cyclic exchange is shown to produce such a confinement. Dopants 
on different chains experience an effective space-extended non-frustrating pairwise spin interaction. This effective 
interaction between impurity-spins is long-ranged and therefore is expected to play a crucial role in the mecanism 
of antiferromagnetic (AF) long-range ordering (LRO) observed in spin-Peierls (SP) compounds such as CuGeOs 
doped with non-magnetic impurities. 
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Fig. 1 . Scliematic picture of the coupled chains model with 
nearest neighbor (NN), next-nearest neighbor(NNN), in- 
ter-chain and 4-spin couplings Ji, J2 = oJi, Jx, and 
J4. Full (resp. open) circles stand for spin-^ sites (resp. 
non-magnetic dopants). 

1. Introduction 



Doping a spin liquid system with non-magnetic 
impurities leads to very surprising new features. 
For example in the doped quasi one-dimensional 
compound Gni-JArcGeOz (M=Zn or Mg), the 
discovery of coexistence between dimerization and 
AF LRO at small impurity concentration x has 
motivated extented experimental [1] and theo- 
retical [2,3,4,5,6,7] investigations. The impurity- 



induced AF LRO has been observed in other 
doped spin liquid materials such as the 2-legs lad- 
der Sr(Cui_a:Zna;)203 [8], the Haldane compound 
Pb(Nii_a;Mga;)2V208 [9] Or the coupled spin dimer 
system TlCui_,Mg,Cl3 [10]. 

Replacing a spin-^ in a spontaneously dimerised 
(isolated) spin chain by a non magnetic dopant 
(described as an inert site) liberates a free spin 
i, named a soliton, which does not bind to the 
dopant [2] . On the other hand, a static bond dimeri- 
sation produces an attractive potential between 
the soliton and the dopant [2,3] and consequently 
leads, under doping, to the formation of local mag- 
netic moments [2,5] as well as a rapid suppression 
of the spin gap [4]. However, a coupling to a purely 
one-dimensional (ID) adiabatic lattice [6] does not 
produce confinement in contrast to more realistic 
models including an elastic inter-chain coupling (to 
mimic 2D or 3D lattices) [6,7]. 

Frustration and inter-chains effects are nec- 
essary to understand the impurity-induced AF 
ordering in the doped spin-Peierls (SP) material 
Cui_2:M2.Ge03. In section 2 we report numer- 
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ical studies of models for doped coupled spin 

chains [11] and concentrate on the local moment 
formation induced by the doping. Dopants on 
different chains experience an effective space- 
extended non-frustrating pairwise spin interaction 
which is long-ranged and therefore is expected to 
play a crucial role in the mecanism of AF LRO. 
In section 3, we report exact diagonalisation (ED) 
results for the effective magnetic coupling which 
appears between released spins. This long-distance 
interaction between impurity-spins is finally in- 
cluded in an effective 2D model with a small con- 
centration X of spins- i put at random on a square 
lattice. Concluding remarks are given in Section 4 
where we also mention some preliminary results 
obtained by a Quantum Monte Carlo (QMC) 
study of an effective diluted model. 



2. Impurity induced local moment 
formation in doped coupled frustrated spin 
chains 

Let us first consider a model of coupled frus- 
trated spin-i antiferromagnetic chains (see Fig.l). 
Following Schulz [12], a mean-field (MF) treat- 
ment of the inter-chain couplings has been per- 
formed [11] and the resulting Hamiltonian is given 

by 

+aSi^a ■ Si+2,a + K,aSla + Constant , (1) 

where 

Ka = Jd{Sl,+ ,) + {Sl,-,)) (2) 

accounts for first-order effects in the inter-chain 
magnetic coupling Jj_, and 

^Ji,a = J4,{{Si,a+l-Si+l^a+l) + {Si,a-l-Si+l^a-l)} 

(3) 

takes a generic form because it might have mul- 
tiple origins; although a four-spin cyclic ex- 
change [13,14] provides the most straightforward 
derivation of it [11], J4 can also mimic higher or- 
der effects in J± [15] or the coupling to a 2D (or 
3D) lattice [6] . i is a lattice index along the chain 
of size L and a labels the M chains {L and M cho- 
sen to be even) . Periodic boundary conditions will 
be assumed in both directions. The energy scale is 
set by the coupling along the chains J = 1 and a 
is the relative magnitude of the NNN frustrating 
coupling. 

In the pure case (i.e. without impurity), all the 
chains are equivalent and the problem is therefore 
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Fig. 2. SP-AF phase diagram in the (a, Jj_) plane from ED 
of chains of length up to 16 sites. Symbols eorrespond to 
different values of J4 > as indicated on plot. Typically, 
FSE are smaller than the size of the symbols. The computed 
transition lines are extended by tentative transition lines 
(dashed lines) in the region where FSE become large. At 
J4 = we have plotted a few points in the vicinity of the 
MG point [16] for L = 12 and L = 16 (reprinted from 
Ref.[ll]). 

reduced to a single chain problem in a staggered 
magnetic field hi = —2J± < > and with its 
NN exchange modulated by 6Ji — Ji{Si ■ Si+i) 
if J4 < or SJi = JiiSi+i ■ Si+2) if Ji > 0. 
Using Lanczos ED up to the convergence of the MF 
procedure [17], we can identify two different phases 
in the {a,J±) plane. A dimerised SP phase and an 
AF ordered phase separated by a transition line 
Jj_ = J±{a) (sec Fig. 2). ED have been performed 
on small systems {L < 16) for different values of J4 . 
Fortunately, the finite size effects (FSE) are small 
in the gapped regime. The modulation created by 
J4 stabilizes the SP phase, as we can observe on 
Fig.2 

Let us now turn to the doped case. A non- 
magnetic dopant is described here as an inert 
site decoupled from its neighbors. Under doping 
the system becomes non-homogeneous so that we 
define a local mean staggered magnetization, 

Mf:' = \i-iy+"{2{SlJ - {SU.^J - {S!_,J). 

(4) 

Following the method used in Ref. [7], the MF 
equations are solved self-consistently on finite L x 
M clusters and lead to a non-uniform solution. At 
each step of the MF iteration procedure, we use 
Lanczos ED techniques to treat exactly (although 
independently) the M non- equivalent timte chains 
and compute (iSf^) for the next iteration step un- 
til the convergence is eventually achieved. We first 
consider the erase of a single dopant. Whereas in 
the case J4 = the soliton remains de-confined 
as can be seen from Fig. 3, a very small J4 7^ 
is sufficient to produce a confining string which 
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Fig. 3. Local magnetization for L X Af =16 X 8 cou- 

pled chains with one dopant D (shown by arrow) located at 
a = 1, i = 16 in the dimerised phase (a = 0.5, J±_ = 0.1). 
Circles correspond to J4 = (shown up to the third 
neighbor chain of the doped one) and squares (crosses) to 
J4 = 0.01 (J4 = 0.08). The coupling J2 across the dopant 
has been set to for convenience (reprinted from Ref. [11]). 

binds the soliton to the dopant. Note that the inter- 
chain couphng induces a " polarization cloud" with 
strong antiferromagnetic correlations in the neigh- 
bor chains of the doped one; we can therefore de- 
fine a typical length scale in the transverse direc- 
tion S^±_ which is ~ 1 in the case Jj^ = 0.1, as we 
will study in the last part of next section. A con- 
finement length in the chain direction ^|| can also 
be extracted. Defined by 



(5) 



we have calculated it for a 16 x 8 system with a — 
0.5 and Jj_ = 0.1, and we show its variation as a 
function of J4 in Fig. 4. FSE decrease for increasing 
J4. Note that i\\{JA) 7^ C|| (~>^4) and a power law [3] 
with different exponents rj is expected when J4 
0. A fit gives r/ - 0.33 if J4 < and r/ ~ 0.50 for 
J4 > (Fig. 4). This asymmetry can be understood 
from opposite renormalisations of Ji for different 
signs of J4. Indeed, if J4 < then SJi^a > and 
the nearest neighbor MF exchange becomes larger 
than the bare one. Opposite effects are induced by 
J4 > 0. 



3. Effective interaction betv^reen 
impurity-spins 

We now turn to the investigation of the effec- 
tive interaction between dopants. Each impurity 
releases an effective spin ^ , localized at a distance 
^ ^11 from it due to the confining potcnticl set 
by J4. We define an effective pairwise interaction 
J'^^ as the energy difference of the S = 1 and the 
5 = GS. When J''^ ^ E{S = 1) - E{S = 0) is 




Fig. 4. ED data of the soliton average position vs J4 calcu- 
lated for a = 0.5 and J± = 0.1. Different symbols are used 
for L X M = 12 X 6 and 16 X 8 clusters. The long-dashed 
line is a power-law fit (see text). Inset shows the magneti- 
zation profile in the doped (a = 1) chain at J4 = 0.08, ie 
^11 ~ 2.5 (reprinted from Ref.[ll]). 

positive (negative) the spin interaction is AF (fer- 
romagnetic). Let us first consider the case of two 
dopants in the same chain, (i) When the two va- 
cancies are on the same sub-lattice the moments 
experience a very small ferromagnetic J'^^ < as 
seen in Fig. 5 with Aa = so that the two effective 
spins i are almost free, (ii) When the two vacan- 
cies sit on different sub-lattices, Ai is odd and the 
effective coupling is AF with a magnitude close to 
the singlet-triplet gap. Fig. 5 with Aa — shows 
that the decay of J°*^ with distance is in fact very 
slow for such a configuration. The behavior of the 
pairwise interaction of two dopants located on dif- 
ferent chains (Aa — 1, 2,3) is shown in Fig. 5 for 
Aa = 1,2,3 for J4 > 0. When dopants are on oppo- 
site sub-lattices the effective interaction is antifer- 
romagnetic. At small dopant separation J'^^{Ai) 
increases with the dopant separation as the overlap 
between the two AF clouds increases until Ai ~ 
2^11 . For larger separation, J'^^{Ai) decays rapidly. 
If dopants are on the same sub-lattice, solitons are 
located on the same side of the dopants [18] and the 
effective exchange J'^^{Ai) is ferromagnetic and 
decays rapidly to become negligible when Ai > 
2^. The key feature here is the fact that the effec- 
tive pairwise interaction is not frustrating (because 
of its sign alternation with distance) although the 
frustration is present in the microscopic underly- 
ing model. AF ordering is then expected (at T — 0) 
as seen for a related system of coupled Spin-Peierls 
chains [7]. 

Our next step is to fit the numerical data in order 
to derive an analytic expression for Jeff and a long 
ranged non frustrationg effective model for diluted 
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Fig. 5. Magnitude of the effective magnetic coupling be- 
tween two impurities located either on the same chain (a) 
or on different ones (b-c-d) vs the dopant separation Ai in 
a system of size L X M = 16 X 8 with a = 0.5, J± = 0.1, 
and J4 = 0.08. Closed (resp. open) symbols correspond to 
AF (F) interactions. Full lines are fits (see text) 

spins- i on a Lj. X Ly square lattice, 

^^'^ = ^ eri er. Jeff (r 1 - r2 ) S,^ ■ Sr^, (6) 

ri,r2 

with Er = 1 (0) with probability x {1 — x), x being 
the dopant concentration. Using only five param- 
eters, two energy scales and three length scales, 
we can fit ED data with very simple mathematical 
expressions. When Aa = (same chain), J"^ ap- 
proximately fulfills J°^{Ai, 0) = Jo(l - Ai/Cy) for 
Ai even and Ai < ^[j*, and J''^{Ai,0) = other- 
wise. For dopants located on different chains and 
on the same sub-lattice (Ai -I- Aa even) one has, 

J^ff(A^, Aa) = -J, exp(-^) exp(-^), (7) 

while if the dopants are on opposite sub-lattices, 
one gets 

J<=«^(Az,Aa) = j;^exp(-^) (8) 

for Ai < 2(^11 and 

„ , Az — 2f 11 Aa 
r''{Ai,Aa) = -Joexp( ^)exp(--^), 

(9) 

for Ai > 2^11 . The fitting parameters are Jq = 0.52, 
Jq = 0.3, = 17.33, = 2.5 and = 1 in the 
case considered here : a — 0.5, Jj_ = 0.1 and J4 = 
0.08 (see Fig.5). 



4. Conclusion 

We can conclude this study by mcntionning 
some preliminary results obtained by the way of 
QMC simulations [19] performed on the effective 
diluted model Eq.(6) with a great number of spins 
Ns < 256 . Even at very small concentrations x, 



a Neel type AF LRO at T = is observed as a 
result of the simulations; details about this study 
will be reported elsewhere [20] . 
We gratefully acknowledge Anders W. Sandvik for 
the interest he took in this work. 
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